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Abstract: 

Space-borne observations of COz from the Orbiting Carbon Observatory-2 are used to 
characterize the response of the tropical atmospheric CO2 concentrations to the strong El Nifio 
event of 2015-2016. Correlations between atmospheric CO growth rate and the El Nifio 
Southern Oscillation have been well known; however, the magnitude of the correlation and the 
timing of the responses of oceanic and terrestrial carbon cycle remain poorly constrained in 
space and time. Here we use space-based CO> observations to confirm that the tropical Pacific 
Ocean does play an early and important role in modulating the changes in atmospheric COz 
concentrations during El Nifio events — phenomenon inferred but not previously observed due to 


lack of high-density, broad-scale CO2 observations over the Tropics. 


One Sentence Summary: 
NASA’s OCO-2 mission provides a first-hand look at the space-time evolution of tropical 


atmospheric CO> concentrations in response to the 2015-2016 El Nifio 
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Main Text: 

El Nifio Southern Oscillation, or ENSO, is the dominant mode of tropical climate 
variability on interannual to decadal timescales (/-5) and is correlated with large inter-annual 
variability in global atmospheric CO, concentrations (6-/9). Studying the response of the carbon 
cycle to this natural climate phenomenon is critical to understand and quantify the sensitivity of 
the carbon cycle to climate variability, and by extension to climate generally (20). Although the 
ENSO cycle originates in the equatorial Pacific, its impact on the carbon cycle is felt globally 
due to its regional teleconnections (22-23) and influences on atmospheric and ocean circulation, 
precipitation, temperature, and fire emissions (/, 24-26). Partitioning the response of the 
constituent components of the carbon cycle to a complete El Nifio event has been challenging 
because of the limited number of CO> observations over the tropical land and ocean regions. 

Observations of atmospheric CO, from space provide a global view of the carbon cycle 
that can be used to describe phenomena that have been previously pieced together from sparse in 
situ data. NASA’s Orbiting Carbon Observatory-2 (OCO-2) mission was successfully launched 
on July 2, 2014 and started providing science data in early September 2014 (70). Within the first 
two years of operation of the OCO-2 mission, a major El Nifio (warm phase of the ENSO) 
occurred (27-30). We provide an approach for studying the temporal sequence of El Nifio- 
induced changes in global CO» concentrations using observations from the OCO-2 mission that 
are validated with the Tropical Atmosphere Ocean (TAO) mooring CO? data. We see a response 
from the tropical Pacific Ocean during the early stages of an El Nifio event and a lagged (and 
much larger) terrestrial signal as the El Nifio reaches maturity. 


El Nino and the global carbon cycle 


Correlations between the atmospheric CO growth rate and El Nifio activity have been 
reported since the 1970s (6-8, 31-32), although the magnitude and timing of the responses of the 
ocean and terrestrial components remain poorly constrained (33). Here, the word terrestrial 
includes both changes in biospheric productivity (respiration and photosynthesis) as well as 
biomass burning (fires). Following previous strong El Nifio events (for example, the 1982-1983 
and 1997-1998 El Nifio events), methods for measuring the atmospheric CO2 response to ENSO 
were based on in situ atmospheric CO? observations at a handful of surface stations that transect 
the tropical Pacific, including Mauna Loa, Christmas Island and American Samoa (8, 34) as well 
as shipboard transect measurements (/2, 35-36). The annual growth rate of atmospheric CO2 
measured at these remote stations and other sites around the globe show remarkable correlation 
with ENSO indices, with a rapid increase in atmospheric CO2 associated with the late stage of an 
El Nifio event (/9, 37). The ocean response to El Nifio events is based on studies looking at in 
situ observations, for example, surface ocean pCO? observations from ships of opportunity (/2), 
moorings (38-39), or targeted field campaigns during El Nifio events (9-10, 40-41), and a variety 
of mechanistic ocean models (24, 52, 54, 61, 65, 67). 

The overall increase in the release of CO> to the atmosphere during strong El Nifio events 
has been attributed to a decrease in biospheric uptake of COz (e.g., due to drying of tropical land 
regions and an increase in plant and soil respiration) combined with enhanced fire emissions. In 
recent years, this has led to a growing body of literature (42-49) concluding that ENSO-mediated 
variability in tropical net land primary productivity is what primarily influences the atmospheric 
CO» growth rate. A handful of studies (25, 50-57) have disputed any consistent or coherent 
response from the land component during El Nifio events, thus highlighting the high level of 


uncertainty and disagreement within the carbon cycle community. 


The El Nifio-CO? signature should have a significant tropical Pacific Ocean component 
as well, with opposite sign to the terrestrial response (/0, /3, 33). During strong El Nifio events, 
there is a large-scale weakening of the easterly trade winds and suppression of eastern equatorial 
Pacific upwelling (indicated by a deeper thermocline) that reduces the supply of cold, carbon- 
rich waters to the surface (Fig. 1). This reduces the usual strong outgassing of CO. from this 
region (52-67), typically on the order of ~0.4—0.6 PgC yr" to the atmosphere, by ~40-60% 
during an El Nifio event (9-/2, 33, 36, 60, 73). If net fluxes were to remain constant elsewhere, 
these substantial net air-sea CO2 anomalies should lead to a reduction in the growth rate of 
atmospheric COs, at least during the early stages of El Nifio. 

Understanding these variations in atmospheric COz, the timing of these variations and the 
underlying processes that cause them have been of great interest within the carbon cycle 
community (/, 10-13, 15, 20, 33, 50). Integrating information from ocean- and atmosphere-based 
estimates, and modeling studies, we now know that it is the combined and opposite effect of 
ocean and terrestrial responses, which contribute to El Nifio-related variations in atmospheric 
CO? (33). What remain controversial though are the timing of the ocean response and a precise 
quantification of its role. This is of crucial importance because typically the interannual 
variability (IAV) in the growth rate of atmospheric CO> is used to constrain the climate 
sensitivity of land carbon fluxes (11) (20-21); however, if a component of the IAV is being 
modified by ocean fluxes, then these inferences of Yr need to be reconsidered. 

Because of the few surface CO2 monitoring stations over the center of action (i.e., 
tropical Pacific Ocean), it has been challenging to directly observe the timing and changes in flux 
of CO, from the ocean to the atmosphere that affect the atmospheric CO2 growth rate during an 


El Nifio event. Efforts to analyze the data from distant measurement locations tend to identify the 


enhanced CO, fluxes from the terrestrial carbon cycle, which dominate during the later stages of 
El Nifio. The high-density, broad-scale observations of CO2 from OCO-2 provide a valuable tool 
to partition the ocean and terrestrial carbon cycle responses to El Nifio. 

Time series of Xco2 anomalies during the 2015-2016 El Nitto 

OCO-2 observations describe the column-averaged CO, dry air mole fraction (Xco2). 
More details regarding the OCO-2 mission, data features, Xco2 retrievals, etc. are provided in the 
Supplementary Materials, and are available in (69) and (70) while validation of Xco2 via 
comparisons to a ground-based network are provided in (7/). 

El Nifio events are identified by warm sea surface temperature anomalies in precise 
regions of the tropical Pacific Ocean, with the most commonly used being the Nifio 3.4 region 
(5°S-5°N, 170°W-120°W). Figs. 2A and 2B show the trend in Xco2 anomaly (90) for the Nifio 
3.4 region and its temporal evolution relative to two ENSO indices (9/), including the Oceanic 
Nifio Index (ONI - derived from sea surface temperature anomalies in the Nifio 3.4 region) and 
the Southern Oscillation Index (SOI - derived from observed sea level pressure differences 
between Tahiti and Darwin, Australia). The 2015-2016 El Nifio began around March 2015 and 
reached its peak over the Central Pacific between November 2015 and January 2016 (30). The 
Xco2 anomaly (Fig. 2B) shows two distinct periods over the entire El Nifio event: (a) 
Development phase of El Nifio (Spring-Summer 2015) — we argue that the negative Xco2 
anomaly is due to a reduction in local CO2 outgassing from the tropical Pacific Ocean, and (b) 
Mature phase of El Nifio (Fall 2015 onwards) — we argue that the positive trend in Xco2 anomaly 
reflects an increase in atmospheric CO, concentrations due to terrestrial sources (i.e., 
combination of reduced vegetation uptake across pan-tropical regions and enhanced biomass 


burning emissions from SE Asia and Indonesia). The time series in Fig. 2B shows the space- 
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based CO, dataset documenting the response of the carbon cycle (both oceanic and terrestrial) 
during an entire El Nifio event, capturing both the development and the mature phase and the 
transition between those two. The timing of the OCO-2 launch was extremely fortuitous in this 
regard. 

Deriving the Xcoz anomalies require observations taken by both NASA’s OCO-2 and the 
Japan Aerospace Exploration Agency’s (JAXA) Greenhouse Gases Observing Satellite 
(GOSAT) (68) mission. The short OCO-2 record makes it impossible to fit a long-time series and 
calculate anomalies, and hence data from the GOSAT mission (operating since January 2009) 
was utilized to generate the Xco2 climatology. The OCO-2 team retrieved Xco2 from the first 7 
years of the GOSAT observations using the same retrieval algorithm that generated the OCO-2 
data product (90). Continuous global coverage from these two missions allows us to stitch 
together a long-time series of Xco2 over remote regions, such as the tropical Pacific Ocean (Figs. 
S1-S2). However, utilizing two data sources, i.e., GOSAT and OCO-2, can incur errors in the 
analyses due to changes in the two instruments, their observing strategies and sampling density. 
Fig. 2B also illustrates the corresponding uncertainty in our analyses. The uncertainty is 
calculated using an ensemble technique (Section C in Supplementary Materials) and further 
brings out the two phases in the time series of the Nifio 3.4 Xco2 anomaly — 0.3 ppm 
uncertainties during the El Nifio development phase with both the upper and lower bounds below 
the zero line, and larger uncertainties of +0.5 ppm during the mature phase of the El Nifio event. 
These larger uncertainties during the latter stages of the El Nifio illustrate the challenge in 
attributing the changes in Xcoz2 anomalies to the competing, and often opposing, signals from the 


ocean and the terrestrial components of the carbon cycle. 


Attributing the two observed phases of Xco2 anomalies to the ocean and the terrestrial 
response 

Our argument for the two observed phases in the Xco2 anomaly time series is supported 
by complementary data sources. The ocean response is corroborated by sea surface pCO2 
observations from an in situ network of autonomous CO2 systems on the TAO moored buoy 
array (9, 38, 72). Data are not directly comparable to atmospheric Xco2 as they describe CO 
variations at the ocean surface. The trend of the difference between the sea surface and 
atmospheric COz (ApCOz), however, does capture typical El Nifio signatures. For example, Fig. 
2C illustrates data from one of the moored buoys in the Nifio 3.4 region (0°, 170°W), which 
shows decreasing ApCO> over the spring and near-zero ApCOz by December 2015. A 
suppression in the upwelling of CO-rich waters caused by weakening of the easterly trade winds 
leads to a reduction in the surface ocean carbon content, which in turn leads to a decline in the 
magnitude of sea-to-air CO> fluxes. The flux estimates at this buoy location are 1.35 + 0.21 (1e) 
gC m~ month” during the November 2014 to February 2015 period (i.e., non El Nifio 
conditions) that gradually decrease to 0.087 + 0.083 (1c) gC m™ month” between November 
2015 and February 2016 (1.e., El Nifio conditions). This indicates a near-total shutdown of sea- 
to-air flux during Boreal Winter 2015-2016 relative to the neutral 2014-2015 Boreal Winter. 
Previous studies focusing on the tropical Pacific Ocean have reported flux reductions of ~40- 
60% over the entire basin (9-/2, 33, 36, 60, 73). Atmospheric transport model calculations with a 
prescribed set of flux patterns and comparing to the observed Xco2 anomalies (Section A in 
Supplementary Materials) suggest a flux reduction of ~26-54%. 

While these numbers are roughly similar, we do recognize the limitation in comparing 


flux estimates from one point (namely the TAO location at 0°, 170°W) to flux estimates for the 


entire Nifo 3.4 region and/or the tropical Pacific Ocean from previous studies. Large-scale 
changes in the physical and biogeochemical dynamics during El Nifio events result in significant 
spatial and temporal variability in the surface pCO> distributions (/2, 62, 65). Additionally, these 
spatial variations and their seasonal progression are uniquely tied to each El Nifio event; thus, 
different flavors of El Nifio events and/or shifts in the El Nifo phenomena (86-88) will influence 
the evolution of the seasonal cycle of pCO> and air-sea CO? fluxes over the region. For the 2015- 
2016 El Nifio event, the TAO buoy at 0°, 170°W lay closest to the edge of the warm pool and 
registered the first response to the onset of El Nifo conditions. As observations from other TAO 
locations (92) are becoming available, it is evident that in the eastern part of the basin there was 
an overall suppression of the outgassing CO> source but with large variability in pCO». Data 
synthesis and modeling work with these and other in situ observations are ongoing to quantify 
the exact magnitude of ocean CO? fluxes over different tropical Pacific regions during the 2015- 
2016 El Nifio. 

The second phase in the Xco2 anomaly time series is driven by the terrestrial component 
of the carbon cycle, and the transport of this signal to the remote Nifio 3.4 region. The anomalous 
increase in CO can be attributed to a combination of terrestrial sources, including a reduction in 
the global biospheric uptake, increases in soil and plant respiration and enhanced fire emissions. 
In fact, the impact of enhanced fire emissions and their regional progression was a well-studied 
feature following the strong 1997-1998 El Nifio (26, 43, 74-76). For the 2015-2016 El Nifio 
event, strong correspondences between Xco2 from OCO-2 and the carbon monoxide (CO) total 
column anomalies from the Measurements of Pollution in the Troposphere (MOPITT) 
instrument on the NASA Terra platform, are evident over the tropical Pacific Ocean, especially 


during Fall 2015 (Fig. 2D). We conjecture that these CO total column anomalies are 


representative of the emissions from the 2015—2016 Indonesian peat fires (77-80), which were 
advected into the tropical Pacific region. El Nifio-related changes in the Walker circulation (1.e., 
westerly winds) and the slightly more southern than normal positioning of the Inter Tropical 
Convergence Zone (ITCZ) (8/) may allow emissions from the Indonesian peat fires to carry over 
into this region (Fig. S4). It is interesting to note from Figs. 2B and 2D that the positive increase 
in Xco2 anomaly actually leads the fire signals by 1-2 months. This indicates that the release of 
carbon flux resulting in an increase in CO? concentrations is only partially pyrogenic; reduced 
vegetation uptake due to droughts is a significant contributor, and quite possibly the initial cause 
of the increase in Xco2 anomaly. 
Isolating the observed negative Xco2 anomaly to an ocean signal 

The time dependence of the Xco2 anomalies during the 2015-2016 El Nifio indicate that 
the initial decrease in atmospheric CO: is due to suppression of upwelling in the tropical Pacific. 
This early negative response is subsequently offset by a large positive anomaly due to the 
terrestrial component. Assuming no significant interannual changes elsewhere in the global 
ocean, we can further confirm our argument by a comparison of the Xco2 anomaly in the Nifio 
3.4 region with the global Xco2 anomaly (Fig. 4A). By differencing the far-field effect from the 
local signal, the influence of the reduction in CO, outgassing from the tropical Pacific Ocean is 
clearly visible during the onset phase of El Nifio. The peak reduction registered over the Nifio 3.4 
region relative to the global Xco2 anomalies is 0.35 ppm in June 2015, which occurs a couple of 
months after the initiation of the El Nifio event. Lag correlation of the Nifio 3.4 Xco2 anomalies 
against the ONI index indicate that the highest positive correlation occurs when the 
concentration-related anomalies lag the SST-related anomalies by 1-2 months (93) (Fig. S8). The 


time lag relationship can be precisely quantified during the onset phase of El Nifio, but it is much 
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more difficult to interpret during the succeeding El Nifio stages when any reduction in CO, from 
decreased equatorial upwelling is masked by the signal from terrestrial processes. Thus, if it were 
not for the reduction in outgassing from the ocean, the impact from the terrestrial sources would 
likely be larger. Our analysis confirms the findings from (/3) that the slowdown of atmospheric 
CO, increase during the early stages of an El Nifio is indeed related to the decreased sea-to-air 
flux of CO> in the tropical Pacific Ocean. The coverage from the OCO-2 mission has enabled us 
to verify this hypothesis and monitor its temporal evolution using real atmospheric CO? 
observations. 

The early stage negative Xco2 anomaly is unique to the tropical Pacific Ocean and is not 
influenced by global, terrestrial or large-spatial scale fluxes. Due to the large interhemispheric 
gradients in CO, typical variability in tropical CO. concentrations can be an aliasing of 
terrestrial processes occurring at higher latitudes. In order to confirm that the recovered ocean 
signal in the Xco2 anomaly is unique to the tropical Pacific Ocean, we examined three other 
ocean regions - the subtropical North Pacific (20°-30°N, 120°-170°W), the subtropical South 
Pacific (20°-30°S, 120°-170°W) and the tropical Atlantic Ocean (5°N -5°S, 5°-35°W). Fig. 3 
shows the specific regions (aside from Nifio 3.4) that we have analyzed, and each of which assist 
us to reject alternative hypotheses. Non-zero differences in Xco2 anomalies between these and 
the Nifio 3.4 region (Fig. 4) indicate that the trend observed over the tropical Pacific Ocean is 
distinct from other ocean basins. This makes intuitive sense from our mechanistic understanding 
as well - while large impacts of ENSO on the sea-to-air COz flux in the tropical Pacific Ocean 
are expected, studies have shown minute and delayed influence of the ENSO modes on the 
variability of carbon fields in the tropical Atlantic Ocean (66, 82-83). 


Perspective 
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The strong El Nifio in 2015-2016 caused a reduction in the magnitude of CO, outgassing 
from the tropical Pacific Ocean. These changes, albeit of varying magnitude, extended over a 
large portion of the tropical Pacific, and impacted the large-scale modulation of the physical 
processes responsible for the CO efflux from this region. Almost all observing networks (i.e., 
OCO-2, TAO, etc.) were aided by the strength of this signal. However, OCO-2 provided a more 
comprehensive view of the tropical Pacific Ocean signal than previous observing networks given 
its: (a) greater coverage and more frequent sampling than in situ networks, and (b) improved 
resolution and precision than earlier space-based instruments. For example, GOSAT, like OCO-2 
is sensitive to the total CO2 column, but has lower precision (2 ppm single sounding random 
error for GOSAT vs. 0.5 ppm for OCO-2) and lower sampling density (100x less soundings). 
The immediate next step will be to fold in these observations into an inverse modeling 
framework (/3, 15, 50, 56) to infer the underlying net fluxes between the ocean and atmosphere 
and the terrestrial biosphere and the atmosphere. This would help establish the real benefit of 
OCO-2, especially against the backdrop of previous studies that had to rely on sparse 
atmospheric constraint to infer changes in CO> surface fluxes during El Nifio events. 

Based on OCO-2 data alone, however, we cannot quantitatively discriminate the relative 
roles of reduction in biospheric activity uptake due to a warmer and drier climate in 2015 versus 
enhanced fire emissions. While we can quantify the temporal response of the ocean versus the 
terrestrial component and qualitatively observe the gradients in the response of different tropical 
Pacific Ocean regions (Fig. 5), it is much more challenging to discriminate the contribution of 
fire emissions and the delayed response of the terrestrial biosphere to El Nifio-induced changes 
in weather patterns. The impact of ENSO is typically felt by the terrestrial biosphere over several 


months to a year after the actual event. Studies on both progressions of droughts (84) and fires 
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(26) during an El Nifio cycle have shown a hysteresis in the Earth system’s response to changes 
in temperature and precipitation patterns. Analyses using ancillary data sources such as solar- 
induced fluorescence (SIF), bottom-up model simulations and inverse modeling calculations are 
typically necessary to quantify the partitioning of the terrestrial carbon fluxes (reduction in 
biospheric uptake vs. increase in fire emissions) as has been pursued in a companion study (85). 

Our study provides a short-term perspective on the potential of CO2 observations from 
space for unraveling more complex relationships of carbon sources and sinks in the future. A 
longer time series of observations will enable testing more hypotheses such as the possibility of 
regionally dependent gradients in air-sea CO fluxes in the tropical Pacific, or adding data to 
support biogeochemical theories at previously inaccessible scales. From a long-term perspective, 
such information will improve our process-based understanding, inform our current suite of 
mechanistic models, and ultimately, better constrain future carbon cycle projections. 
Concluding remarks 

The strong El Nifio event of 2015-2016 provided us with an opportunity to study how the 
global carbon cycle responds to changes in the physical climate system. With the high-resolution 
(both spatial and temporal) observations available from OCO-2, we are able to directly: (a) 
observe the strong correlations that exist between atmospheric CO2 concentrations and the El 
Nifio signal, and (b) track the development of the atmospheric CO2 anomaly as it switches from a 
negative phase (i.e., due to a reduction in CO> outgassing from the tropical Pacific Ocean) to a 
strong positive phase (i.e., due to a reduction in biospheric uptake and increased fire emissions). 
The most important contribution of the space-based OCO-2 mission is the ability to observe and 
monitor carbon cycle phenomena at high-density over large spatial scales, which has not been 


possible from the existing in situ network. 
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The complexity of the El Nifio — CO? signature illustrates that it is a multifaceted system 
with contributions from many regions and processes. Understanding and predicting its behavior 
requires separating out the many terrestrial and marine regions that contribute (/, 33) and 
identifying both the geophysical (3, 27, 30) and the biological (/0, 59, 89) phenomena that 
respond in their own unique ways. However, the impact on the carbon cycle is unified through 
the global mixing of CO, in the atmosphere - OCO-2 makes a unique contribution by providing 
both the global coverage and fine surface spatial detail; alongside the in situ CO2 network of 
moorings and shipboard measurements provide the long-term climate-quality record of 
atmospheric and ocean CO, observations and serves to validate the OCO-2 observations and 
model products. We emphasize that this diverse observing portfolio is necessary, and the 
complementary information provided by these observing systems will likely prove critical in 
understanding the partitioning of carbon fluxes during the 2015-2016 El Nifio, the relative 
contribution of ocean vs. land to the global atmospheric CO, growth rate, and the sensitivity of 


the carbon cycle to climate forcing on interannual to decadal timescales. 
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Fig. 1. Schematic of the mechanistic differences between normal (A) and El Nifio (B) conditions 
and associated carbon response over the tropical Pacific Ocean. Warm ocean surface 
temperatures are denoted in red and cooler waters in blue. During El Nifio conditions, easterly 
trade winds weaken and westerly wind bursts occur. In association with the shift in wind 
regimes, the western tropical Pacific warm pool moves eastward and the slope of the thermocline 
flattens in the central and eastern tropical Pacific. This suppresses upwelling of cold, carbon-rich 
waters in the central and eastern tropical Pacific, reducing the magnitude of CO, outgassing into 
the atmosphere. Also shown are changes in atmospheric convection, wherein convection shifts 


eastward in response to eastward displacement of western tropical Pacific warm pool waters. 
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Fig. 2. OCO-2 observes the response of the carbon cycle for an entire El Nifio event. Temporal 
evolution of (A) the 2015-2016 El Nifio as captured by the ONI and the SOI indices, (B) Xcoz2 
anomalies and associated uncertainties in the Nifio 3.4 region, (C) ApCOz2 from the TAO 0°, 


170°W mooring, and (D) the CO total column anomalies in the Nifio 3.4 region. 
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Fig. 3. Schematic showing the specific ocean basins (Nifio 3.4, N. Pacific, S. Pacific and Trop. 
Atlantic) that were analyzed in this study. Xco2 anomalies are calculated for these different ocean 
basins, and subsequently compared to the Xco2 anomalies from the Nifio 3.4 region. Each of 
these regions was considered to accept/reject a specific hypothesis that could potentially bias the 
observed trend in the Nifio 3.4 Xco2 anomalies. After rejecting these hypotheses, we conclude 
that the negative Xco2 anomaly observed over the Nifio 3.4 during the onset phase of El Nifio 


2015-2016 is unique and has to be driven by local changes in the ocean fluxes. 
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Fig. 4. Difference in Xco2 anomalies between the Nifio 3.4 region and (A) the globe, (B) the 
tropical Atlantic Ocean, (C) the subtropical Pacific Ocean from September 2014 to May 2016. 
Definitions of the regions are provided in Fig. 3. In Panel (A), we see a robust pattern of negative 
Xco2 anomaly between Nifio 3.4 and the globe that is largest in 2015 and well synchronized with 
the onset phase of El Nifio. In Panels (B) and (C), non-zero differences between Nifo 3.4 and the 
other ocean basins indicate that the Nifio 3.4 trend is not reproducible in other ocean basins; thus, 
allowing us to attribute the negative anomaly in Fig. 2B to a reduction in local CO> outgassing 


over the tropical Pacific Ocean. 
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Fig. 5. Time evolution of the Xco2 anomalies (ppm) averaged over 5°S-5°N. The x-axis 
represents longitude and the y-axis shows the time progressing from top to bottom in months. 
The 2015-2016 El Nifio event and its onset and mature phases are highlighted to show the 
distinct responses observed over the tropical Pacific Ocean. The grey dashed lines capture the 
boundaries of the Nifio 3.4 region. During the onset phase (i.e., March — July 2015), perceptible 
gradients are observable from the far western Pacific to the central Pacific (consistent with the 
increasing flux from west to east) along with high variability in the Xco2 anomalies in the central 
Pacific. We also notice that the Xco2 anomalies are smaller over the eastern Pacific, which is 
consistent with surface seawater pCO? data collected on the TAO buoys (92). The transition 
from the ocean to the terrestrial signal happens between July and October 2015. Towards the 
latter stages of the El Nifio event (i.e., November 2015 and later), the terrestrial signal dominates 


the observed trends in Xcoz likely masking any underlying ocean signal. 
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Materials and Methods 


Xcoz retrievals from OCO-2 and GOSAT-ACOS 


OCO-2 is NASA’s first dedicated satellite mission for measuring column average, 
atmospheric carbon dioxide (CO2) dry air mole fraction (Xcoz2) with the accuracy, resolution, and 
coverage needed for quantifying CO> fluxes (sources and sinks) on regional scales over the globe 
(69, 94-95). Over the sunlit hemisphere, the OCO-2 mission collects nearly one million 
soundings per day (24 soundings/second) at approximately 3km/ nadir resolution. After cloud 
and aerosol screening, between 7 and 12% of these soundings yield full-column estimates of CO2 
on monthly timescales (96-97). 

The OCO-2 observing strategy (nadir/glint modes) has been refined over the first two 
years of operations to improve the measurement coverage and yield, especially over the ocean. 
The initial observing strategy recorded only glint or nadir observations over the entire sunlit 
hemisphere for a complete, 16-day, ground-track repeat cycle, and then used the other observing 
mode in the next 16-day cycle. This approach provided adequate coverage of oceans and 
continents on monthly timescales. In July 2015, this observation strategy was modified to 
alternate between glint and nadir observations on alternate orbits to yield more continuous 
coverage of the entire sunlit hemisphere every day. In November 2015, the observation strategy 
was further optimized to always-collect glint data on orbits that were primarily over the ocean. 

The Xco2 retrievals used in this work are based on the version 7B (v7B) Level 2 
algorithm. These data are freely available via the Goddard Earth Sciences Data and Information 
Services Center (GES-DISC) from the start of mission operation. OCO-2 results are also being 
cross-calibrated and cross-validated with measurements and data products from the Japanese 


Greenhouse gases Observing SATellite (GOSAT, nicknamed "Ibuki"), so that these two satellite 
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datasets can be combined to produce a uniform Xco2 climate data record for use by the carbon 
cycle science community. The GOSAT Xco>? retrievals used in this study have been generated by 
version 7.3 (v7.3) of the ACOS algorithm (GOSAT-ACOS, 98-99). Only the high (H) gain 
observations are retained from the GOSAT-ACOS dataset. These data are bias-corrected using 
the same predictors as those used for the OCO-2 v7B dataset. 

The Xco2 estimates have been validated against results from the Total Carbon Column 
Observing Network (TCCON) and other standards to assess their accuracy and correct regional 
scale biases (7/, 100). After bias correction and data screening, the median residual bias between 
OCO-2 and TCCON Xcoz estimates are less than 0.5 ppm while the root-mean-square (RMS) 
differences between the two estimates are typically less than 1.5 ppm (7/). Direct validation of 
the OCO-2 Xco2 estimates against TCCON sites has only revealed biases beyond 2 ppm for 
some months at Wollongong, Australia (34°S), whereas Ascension Island (10°S in the Atlantic 
Ocean) has not revealed any noticeable biases. Since there are no TCCON stations in the tropical 
Pacific Ocean, we must infer the quality of the ocean data from these validation sites in other 
ocean basins. The OCO-2 team has also noticed errors in the data over the Southern Hemisphere 
ocean that have a seasonal dependency. Those data are screened out using an airmass-dependent 
filter (70). Because the evaluation of systematic errors in these satellite retrievals is currently the 
subject of active research, we have carried out additional sensitivity tests (Section C) to ensure 
that these errors do not impact the findings in this study. 

Typically, the Xco2 retrievals are only performed on scenes nearly devoid of cloud and/or 
optically thick aerosol. Because of this, ENSO-induced changes in cloud patterns may cause 
slight changes in the locations of OCO-2 measurements, but are not expected to induce 


additional biases in these measurements. We find that the retrieval algorithm sees enough clear- 
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sky scenes to consistently generate an adequate number of high-quality soundings for the entire 
time period of this study. For a typical month (e.g., August 2015), we find that GOSAT-ACOS 
returns ~10K “good-quality” soundings while OCO-2 returns ~150K “good-quality” soundings 
globally. For OCO-2, 60% of the retrieved soundings are ocean glint observations. Thus, in 
August 2015, the total number of ocean glint soundings is ~96K. Both of these space-based 
missions provide significantly more coverage over the tropical Pacific Ocean relative to the 
Sparse in situ monitoring network (Fig. S1). In fact, it is the continuous coverage from GOSAT 
and OCO-2 that allows us to generate a time series of Xco2 over the tropical Pacific Ocean (Fig. 
S2) — for each month between June 2009 — May 2016 individual good-quality ocean glint 
soundings over the entire Nifio 3.4 region are aggregated to generate a single monthly Xco2 


estimate over Nifio 3.4. 
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Generation of Xco2 anomalies 


The time-series of Xco2 (and in general, time-series of atmospheric CO2 concentrations) 
exhibit both a linear trend and a cyclo-stationary component due to the seasonal cycle (Fig. S2). 
To account for the seasonality and the upward trend of CO2, we have adopted a two-step 
approach: (1) for each month, all Xco2 retrievals from GOSAT-ACOS and OCO-2 are averaged 
over pre-specified domains (e.g., Nifio 3.4, tropical Pacific Ocean, tropical Atlantic Ocean, 
global) assuming no temporal correlation between months, and (2) for an individual month, say 
August, we find a linear trend that best fits the Xco2 data from 7-years of GOSAT-ACOS and 
OCO-2 observational record for that month, i.e., find the linear trend for August 2009, August 
2010, August 2011, August 2012, August 2013, August 2014, August 2015. Residuals from this 
linear trend (Fig. S3) are the Xco2 anomalies that will be analyzed. Note that the first derivative 
of the trend (slope of the linear regression line) provides an estimate of the monthly Xco2 growth 
rate. Simply comparing the growth rate values with those reported in previous studies delivers a 
necessary (but not sufficient) sanity check of our methodological framework. We have also 
explored the sensitivity of our calculated anomalies with respect to existing fitting methods 
(Section B, Fig. S5) and did not see any significant impact on our findings. 

Mathematically, if €mo, represents the monthly Xco2 anomalies for a given year, then 


these are derived from Equation S1: 


Cmoyr — Vmoyr — Imo,yr = Ymo,yr — (Sina iF Base) nee Equation S1 


where, Ymo,yr= observed monthly Xcoz value for a given year, x 1,..., 7 for the 7 years, and 


ye 
a and f are the intercept and slope parameters that are estimated during the linear regression 


procedure. From a physical standpoint, B». represents the monthly Xco2 growth rate (ppm yr’) 


and @mo represents the monthly Xco2 offset (ppm) relative to the start of our analyses period in 
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2009. This approach provides the flexibility of calculating the linear trend over any spatial 
domain — for example, globally, or over the entire tropical Pacific Ocean or separately for 
individual Nifio regions. On monthly timescales, we do find slight differences in growth rates 
between the entire tropical Pacific Ocean and individual Nifio regions — typically these 
differences are on the order of 0.01-0.2 ppm yr’. On average the global growth rates vary 
between 1.88-2.37 ppm yr’ depending on the season. These numbers are consistent with 
previous studies (/0/-/02) that have attempted to derive space-based Xco2 growth rate over 
various spatial domains (e.g., global vs. Northern Hemisphere vs. Southern Hemisphere). 

Due to the short time series of the OCO-2 record, we had to rely on the GOSAT-ACOS 
dataset to provide a climatological reference. While care was taken to use Xco2 that are based on 
the same retrieval algorithm (OCO-2 v7B and GOSAT-ACOS v7.3), we recognize that stitching 
together data from two different instruments, i.e., GOSAT- and OCO-2 involves an implicit 
change in sampling strategy plus changes in observational density. Within the atmospheric trace 
gas remote-sensing community, this is a strategy that has been employed for generating 
climatologies from data available over a limited timespan — a typical example being the 
generation of ozone climatology using profile data from ozonesondes, observations from 
multiple satellite instruments, etc. (103-105). Section C discusses the sensitivity tests we have 
carried out to confirm that the derived Xco2 anomalies are not being impacted due to this 
‘change of instrument’ issue. In the future, as the OCO-2 data record grows, we will able to 
revise the climatology based on a homogenous Xco2 record, either from actual OCO-2 data or 
quasi-operational Level-3 maps (/06) that will be made available as part of the OCO-2 product 


suite. 
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pCO) observations from the TAO array 


The TAO (Tropical Atmosphere Ocean) array of moored buoys in the tropical Pacific 
Ocean provides real-time, in situ meteorological and oceanographic measurements (72). 
Atmospheric and surface seawater partial pressure of CO2 (pCOz) is currently measured by 
Moored Autonomous pCO (MAPCO>2) systems maintained on the TAO array at 0°, 110°W; 0°, 
125°W; 0°, 140°W; 0°, 155°W; 0°, 170°W; 0°, 165°E; and 8°S, 165°E (38). 

In brief, the MAPCOz system (39) utilizes an automated equilibrator-based gas collection 
system to measure surface seawater xCO> (the mole fraction of COz in air in equilibrium with 
surface seawater) at approximately 14 cm depth and atmospheric xCO> at approximately 1.5 m 
above the sea surface every 3 hours. This xCO2 (wet) measurement is made by a non-dispersive 
infrared gas analyzer (LI-820, LI-COR) calibrated before, during, and after field deployment 
with reference gases traceable to World Meteorological Organization standards prepared by 
NOAA’s Earth System Research Laboratory. The MAPCO> system also measures sample 
temperature, pressure, and relative humidity to calculate xCOz (dry) based on the equations in 
(107). Sea surface temperature (SST) and salinity from TAO temperature and conductivity 
sensors are then used to calculate pCO consistent with ocean carbon standard operating 
procedures as described in (39). 

Once data are recovered from the field and quality controlled, 3-hourly CO, observations 
are archived at the Carbon Dioxide Information Analysis Center 
(http://cdiac.ornl.gov/oceans/Moorings/) and the National Centers for Environmental 
Information (https://www.ncei.noaa.gov/). Monthly averaged real-time observations of ApCO> 
(seawater pCO — atmospheric pCO2) from September 2014 — May 2016 are shown in Fig. 2C 


and Table S1. Quality control checks on the pCO? observations reveal good measurements were 
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collected during this time period. Based on laboratory tests at NOAA Pacific Marine 
Environment Laboratory (PMEL) and field intercomparisons, estimates of uncertainty for 
quality-controlled air and seawater pCO» measurements are <2 atm (39); however, uncertainty 
of preliminary monthly-averaged data are likely slightly higher (<5 atm). Data plots from all 
TAO pCO, locations, which include both real time and finalized data, are available at - 
www.pmel.noaa.gov/co2/story/Opent+Ocean+Moorings. 

Figure 2C captures the temporal trend in the ApCO2 values from the 0°, 170°W TAO 
array buoy, which is located in the western-most portion of the Nifio 3.4 region. Preliminary 
analyses seem to suggest that this buoy registered the first response to the onset of El Nifio 
conditions. On the other hand, observations at 0°, 110°W from the TAO array buoy 
(https://www.pmel.noaa.gov/co2/story/Open+Oceant+Moorings) seem to indicate that upwelling 
of high dissolved inorganic carbon (DIC) water continued in the central and eastern-most portion 
of the Nifio 3.4 region till early Fall 2015 with only brief periods (November 2015 — February 
2016) where surface seawater pCO> was close to equilibrium with atmospheric values. 
Observations from these more eastern sites indicate diverse regimes operating across the tropical 
Pacific; Fig. 5 qualitatively demonstrates that the Xco2 anomalies capture these regional 
gradients depending on its magnitude, seasonality, location, etc. As more in situ measurements 
become available, quantitative investigations are ongoing to check how well the OCO-2 


observations can resolve such regional gradients. 
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CO observations from the MOPITT instrument 


Since March 2000, the MOPITT instrument on board the NASA/EOS Terra platform has 
been monitoring the CO content in the troposphere. Details about the retrieval algorithm, its 
validation and uncertainties are provided in Deeter et al. (108, 112-113) and Emmons et al. (109) 
while analyses and discussion of decadal CO trends are provided in more recent studies (//0- 
111). An instrument cooler failure between May and August 2001 significantly impacted the 
retrieval mean levels before and after the instrument anomaly. Hence for the sake of 
homogeneity, we select data during the period June 2002 - May 2016 for generating a 
climatological value of CO content in the atmosphere. Unlike the shorter Xco2 record, the 
availability of a long and homogeneous CO data record (>14 years) from a single instrument 
makes the generation of the climatology (and associated anomaly calculations) straightforward. 
Based on the recommendation of the MOPITT team, we use the Level 3 MOPITTv6 CO (//2) 
estimated from the thermal-infrared (TIR) channel. TIR measurements offer the best description 
of CO emissions from fires due to the large thermal contrast between the lower troposphere and 
free troposphere due to intense surface heating. For this study, we look at the CO Volume 
Mixing Ratio (VMR) for both the total column and at an individual atmospheric pressure level at 
700 hPa. 

Given that the two largest surface sources of CO are the combustion of fossil fuel and the 
combustion of biomass (forest and savanna fires, biofuel use, and waste burning), MOPITT CO 
anomalies have been used to track emission from fires, including El Nifio-related ones (//4- 
117). While the surface level anomalies can detect the onset of fires immediately (1.e., without 
any time lag), CO anomalies at higher levels in the atmosphere increase gradually and exhibit 


significantly smaller peak anomaly values than the surface. These features are consistent with the 
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expected lag and dilution associated with vertical mixing. Fig. S4 shows that the impact of fire 
emissions between August and October 2015 from S.E. Asia/Indonesia are clearly visible over 
the tropical Pacific Ocean. The largest anomalous sources of CO occur over the Indonesian 
regions, which are emitted as strong pulses over the 3-month window. These maps also rule out 
biomass burning emissions from other regions, such as Africa and South America that can 
potentially influence the rise in CO and CO> concentrations over our region of interest (i.e., Nifio 
3.4). 

Because the signal from enhanced fire emissions are correlated with drought, however, 
we cannot quantitatively discriminate the relative roles of reduction in biospheric uptake due to 
warmer, drier climates or emissions from biomass burning by just looking at the trends in the 
anomalies in the CO2 and the CO data. What we can discern is the timing of the large tropical 
source (of which the biomass burning emissions are a major contributor) and its impact on the 


rise in CO concentration anomalies from August 2015 onwards. 
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Supporting Text 


Section A: Calculation of sea-to-air fluxes from the TAO and the OCO-2 datasets 


Calculation of sea-to-air flux from pCO), observations 

Following Wanninkhof (//8), the net flux of COz across the air-sea interface, F (mass 
area’' time’) is calculated as the product of the gas transfer velocity k (length time’), the 
solubility of CO2, s (mass volume pressure '), and the pCO: (pressure) difference between the 


atm 


ocean and atmosphere (pCO2*” — pCO"), according to Equations S2 and S3: 


F = ks(pC0§” — pC03'™) ... Equation S2 
k = 0.251U2,(S,/660)- _.. Equation $3 


where, Sc is the Schmidt number for CO2, and Ujo is the wind speed at 10m height. By 
convention, the units of k are in cm h”', Uj is in ms" and units of the coefficient 0.251 are (cm 
h'') (ms"')*. Daily averaged wind speeds were acquired from the Modern-Era Retrospective 
analysis for Research and Applications, Version 2, or MERRA-2 (//9). 

Based on these equations, the calculated fluxes at TAO 0°, 170°W are 0.21 + 0.20 (1o) 
gC m” month" between April 2015 — March 2016 (i.e., during the El Nifio 2015-2016 event). 
This is a flux reduction of 84% relative to the neutral April 2014 — March 2015 period during 
which the magnitude of flux was 1.35 + 0.20 (1a) gC m” month”. Note that this flux difference 
is only representative of the western-most portion of the Nifio 3.4 region. The COz fluxes and the 
magnitude of flux reduction could be considerably different in the eastern-most portions of the 
tropical Pacific Ocean. Studies reporting on the magnitude of air-sea flux reduction across the 


tropical Pacific Ocean during El Nifio events typically report a flux reduction of ~40-60%. 
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Calculation of sea-to-air flux from Xco2 observations 

Following Gurney et al. (120), we have used simple atmospheric transport model 
simulations to assess whether the magnitude of ocean flux anomaly based on the pCO2 
observations is consistent with the anomaly observed in the Xco2 estimates. For the 16 models 
participating in the TransCom3 Level 1 experiment (T3L1) (/2/), we have calculated the 
average column mean CO enhancement over the Nifio 3.4 region for the Eastern Equatorial 
Pacific basis function (region 14) fluxes. As per the TransCom experiment protocol, the fluxes 
for this basis function was set to a constant 1 PgC yr’! evenly distributed within the region. These 
data are precomputed and available as part of the TransCom 3 experiments at - 
http://transcom.project.asu.edu/transcom03_output.php. The average difference between Nifio 
3.4 region and a global reference column mean is 0.43 + 0.09 (10) ppm. The global reference is 
based on a latitudinal cutoff of 60°S-60°N keeping in mind the range of latitudes over which 
OCO-2 observations are available. Since the basis function used in these experiments was 1 PgC 
yr’ constant emission, this gives a sensitivity of 1/0.43 = 2.33 PgC yr’ per ppm, or for a 0.26 + 
0.09 (1c) ppm dip (as seen from Fig. 4A over May, June, July, August 2015) about 0.60 + 0.21 
(1c) PgC yr’. Given that the dip in Xco2 anomalies are observed over a period 4 months long 
(Fig. 4A), this would suggest a total flux of 0.20 + 0.07 (1c) PgC. As per the TransCom 
framework, note that this is an implied flux difference, in PgC absolute. Thus, if the original 
magnitude of sea-to-air flux is 0.5 PgC yr (36), then the implied reduction in the flux 
outgassing is 0.2/0.5 * 100 = 40 + 14 (lo) % (1.e., 26-54%). This simple but effective calculation 
provides a necessary sanity check on the magnitude of the signal observed from the Xcoz2 data. 

In reality, the changes in the surface fluxes are far more complicated with processes 


occurring at distinct spatiotemporal scales and with different levels of maturity. The simple 


Dil 


analysis presented here does not take into account spatial gradients in fluxes within the east 
equatorial Pacific region itself or between the east and west tropical Pacific Ocean, all of which 
may have been possible for the 2015-2016 El Nifio event. More detailed analyses using 
atmospheric inverse modeling/data assimilation techniques and both satellite and in situ data are 
necessary (and ongoing) to parse out the grid-scale air-sea flux magnitudes over the tropical 


Pacific Ocean. 
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Section B: Analyzing the sensitivity of Xco2 anomalies to the curve fitting procedure 


Pickers and Manning (/22) discuss that significant bias and uncertainty can be introduced 
in the application of curve fitting programs to atmospheric time series. The specific curve fitting 
method used in this study was designed keeping in mind the attributes of atmospheric CO> 
concentrations (i.e., both long-term growth rate and seasonal cycle). In principle, the overall 
framework is analogous to the one currently used at NOAA (‘CCGCRV)’) to fit observation time 
series from the NOAA in situ network. Originally based on Thoning et al. (123), the CCGCRV 
curve fitting method represents the long-term growth rate by a polynomial function, the annual 
oscillation by harmonics of a yearly cycle and then uses a low pass filter to retain interannual and 
short term variations in the fitted curve. There is no set rule on the number of polynomial terms 
in the function, and it depends on the application and user preference. 

We have used the CCGCRV curve fitting method to detrend our Xcoz time series, and re- 
calculated the anomalies and the growth rates (Fig. $5). A simple linear term was used for the 
polynomial part of the function to make the overall function as comparable as possible to our 
curve fitting method. We find that the overall pattern in the Xco2 anomalies (negative anomaly 
during Summer 2015 followed by a positive anomaly in Fall 2015) is consistent for both 
methods. This implies that the curve fitting procedure did not introduce any spurious bias or 


trend onto our analyses. 
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Section C: Uncertainty analyses for the trend in Xco2 anomalies 


An ensemble approach is used to quantify the error associated with the time series of the 
Nifio 3.4 Xco2 anomaly. Fig. S6 shows specific tracks that GOSAT flies over, and only 6 of 
those tracks lie within the Nifio 3.4 region. We create an ensemble of small regions that are 
roughly 200 km x 400 km around those individual GOSAT tracks. For each of these regions (30 
in total), we apply the same methodology to generate the Xco2 climatology and extract the 
associated anomaly for our study time period (September 2014 — May 2016). The resultant 
spread among the 30 ensemble members defines the lo standard deviation for the mean Xco2 
anomaly line (Fig. 2B), which is still derived from the full Nifio 3.4 region. 

We have carried out an additional sensitivity test (Fig. S7) where we have systematically 
added/subtracted constant mean bias values (0.1-1ppm) to the ocean glint retrievals used in this 
study. Given the lack of validation data over the tropical Pacific Ocean, it is necessary to 
ascertain that residual systematic errors in the ocean glint retrievals (7/, 100) do not alias onto 
our calculations. Fig. S7 illustrates that even in the worst-case scenario, where the ocean glint 
retrievals have a low bias (i.e., negative) of 1 ppm, the overall temporal structure of the Xco2 
anomaly is conserved. The impact on the Xco2 anomaly line is that it shifts up or down 
depending on the magnitude of the bias, which in turn impacts the inference of the magnitude of 
the air-sea fluxes. Certainly, a more complicated structure of the residual biases (e.g., coherent 
spatial and temporal dependency or correlations between biases) will impact these trends in a 
different manner. But based on our current understanding of the satellite retrievals, these 
sensitivity tests and the ocean basin-related diagnostics (Figs. 3 and 4) show that the findings are 
robust to a couple of first order caveats that have not been resolved for the global OCO-2 dataset 


as a whole. 
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Section D: Cross-correlation between Xco2 anomalies and SST anomalies 


We have investigated the cross-correlation between the Xco2 anomalies and the ONI 
index (which are derived from the Nifio 3.4 SST anomalies) for the entire period of June 2009 — 
May 2016 (n=84) to test the magnitude and phase of the relationship between the two variables. 
The largest correlations are found for Xco2 anomalies lagging the ONI by 1-2 months (Fig. S8). 
These correlations are significant at the 95% level based on the random phase test of Ebisuzaki 
(124). Previous studies (6, 13, 14, 31) have looked at carbon cycle related anomalies (i.e., 
anomalies in CO2 concentrations or CO> fluxes) and tested its phasing and strength with a suite 
of ENSO indices (e.g., SOI, Nifio 3 SST anomalies, etc.). Jones et al. (14) provides the closest 
analogy to this work, where they reported a lag of 3 months between the CO2 concentrations 
measured at Mauna Loa and Nifio 3 (5°S-5°N, 90°W-150°W) SST anomalies. 

We believe that the longer response time of the Mauna Loa CO2 anomalies to SST 
anomalies is not necessarily a feature of the carbon cycle but rather due to a time lag associated 
between the driver (SST) and the response (CO2) measured at a location far away from the center 
of action. The availability of Xco2 retrievals directly over the region of action (tropical Pacific 
Ocean and the Nifio 3.4 region) provides a new tool for understanding both the strength of 
correlation of atmospheric CO, changes with ENSO and the real-time lag associated with the 


ocean carbon cycle response. 


oD 


Supporting Figures and Tables 


Ee, * 
0° N RYN 
MARY i \ 
. Okey Y ia vm ENC =m 
RW \\" AN : i Nt nA “| 
rae iN i Ny WE We is WAN \ 
30° S ak ‘ we WA Wi we 
s Nifio 3.4. 


Fig. $1. Coverage over the Pacific Ocean from GOSAT-ACOS v7.3 (magenta) and OCO-2 v7B 
(gold) for one representative month in 2015. Relative to GOSAT, OCO-2 provides more 
observations over the open ocean (~100x more soundings per day) with higher accuracy and 
precision. Also shown are the Nifio regions over the tropical Pacific Ocean that are traditionally 


used to study and define El Nifio events. 
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Fig. S2. Time-series of Xco2 over the Nifo 3.4 region for seven years (June 2009 — May 2016) 


after combining the GOSAT and the OCO-2 data streams. The availability of space-based 


observations over remote locations, such as the Nifio 3.4 region, allows us to stitch together a 


long-time series of Xco2 concentrations from which the Xco2 anomalies can be deduced (Fig. 


S3). Grey-hatched portions indicate time-periods during which data is not available over the 


Nifio 3.4 region (July 2011, June 2012 — July 2012 and June 2014) either due to instrument 


downtime or planned maintenance. 
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Fig. S3. Multi-year time series of Xcoz2 anomalies over the Nifio 3.4 region. The timespan covers 
June 2009 to May 2016. Grey-hatched portions indicate time-periods during which data is not 
available over the Nifio 3.4 region (July 2011, June 2012 — July 2012 and June 2014)- either due 
to instrument downtime or planned maintenance. The horizontal brown lines indicate the range 
of uncertainty (+ 0.5 ppm) expected in the Xco2 data. The vertical dashed red line shows the time 
period from which a subset of this time series (September 2014 to May 2016) has been explored 
in detail in this study. Comparison to previous years illustrates, however, that the anomaly 


signature observed during the 2015-2016 period are distinct relative to other periods. 
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Fig. S4. Spatial distribution of CO total column anomalies retrieved from the MOPITT 
instrument between August and October 2015. The red rectangle shows the location of SE Asia 
and Indonesia. Note the spike in CO anomalies during September 2015 due to the Indonesian 
peat fires, and subsequent advection (visible as a yellow blob) eastward across the tropical 
Pacific Ocean and westward over the tropical Indian Ocean. The peak values of CO anomalies 


occur during October 2015, which is also evident from the time series in Fig. 2D. 
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Fig. S5. Time series of Xco2 anomalies generated by two different curve-fitting methods — the 


procedure used in this study (blue line with diamonds) and the traditional ‘CCGCRV’ procedure 


(orange line with stars). Both methods generate similar trends in the anomaly during the onset 
phase (i.e., negative anomaly in Xcoz) and the mature phase of El Nifio (i.e., steep increase in 
Xcoz anomaly). However, we do note a systematic difference in the magnitude of the positive 
Xco2 anomaly (August — December 2015) likely due to the difference in the number of 


parameters used in the two curve-fitting procedures. 
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Fig. S6. Schematic of the methodology used to generate an ensemble of boxes within the Nifio 
3.4 domain. Each individual box is centered around a GOSAT track — this ensures that the spatial 
coverage provided by the two remote sensing missions (GOSAT and OCO-2) are roughly 
equivalent within the box, thus minimizing the impact of differences in sampling density 


between the two instruments. 


61 


onset phase of peak phase of end of 


El Nifio 2015-16 El Nifio 2015-16 E! Nifo in 
-_————— -_———-} May 2016 
1.5 

—~ even with a 1 ppm bias, 

€ 1p the negative Xco2 =| 

2 anomaly is present 

ym O-9F 

2 

e 0 

re) 

oe -O.5 4+ 

Q 4 anomalies calculated for 

x< [ range of bias values (0.1-1ppm) 

—-1.5 ! L ! | | | ! | | ! 


O D F A J A O D F A 
Months (2014-2016) 


Fig. S7. Uncertainty associated with the Xco2 anomaly time series given residual systematic 
errors in the ocean glint retrievals. A range of constant mean bias values from + 0.1 to + 1.0 ppm 
were added to the ocean glint retrievals from OCO-2, and the Xco2 anomaly re-calculated for 
each bias level. Even in the worst-case scenario (low | ppm bias over the tropical Pacific Ocean), 
the overall temporal structure of the ocean and the terrestrial responses are conserved. The 
magnitude of the Xco2 anomaly shifts upward or downward, which will have consequences for 


the inferred changes in ocean/land —atmosphere fluxes. 
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Fig. S8. Lag correlation of ONI against monthly Xco2 anomalies in the Nifio 3.4 region. Lags 
less than 0 indicate a delay in the response of the Xco2 anomalies to the SST anomalies. The 
horizontal dashed green lines are based on a non-parametric test and indicate the correlation 


required to meet the 95% significance criteria. 
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Month Year Pressure | sw xCO, | air xCQ, SST sss ApCO2z 
(hPa) (umol (umol (°C) (psu) (uatm) 
mol") mol") 
10 2014 1008.4 475.0 397.5 28.42 35.31 74.3 
11 2014 1007.1 471.5 397.1 28.91 35.08 71.1 
12 2014 1005.7 469.3 397.5 29.12 34.72 68.5 
1 2015 1005.9 480.2 398.1 28.68 35.15 78.4 
2 2015 1006.5 483.0 397.9 28.50 35.18 81.3 
3 2015 1006.5 470.0 398.1 28.97 35.15 68.6 
4 2015 1006.7 458.6 397.1 29.59 35.07 58.7 
5 2015 1007.4 431.5 396.7 29.93 34.82 33.2 
6 2015 1006.1 409.8 396.7 30.18 34.54 12.5 
7 2015 1007.8 407.3 396.7 30.02 34.44 10.1 
8 2015 1007.1 402.6 396.5 30.36 34.38 5.8 
9 2015 1007.4 404.9 396.0 30.53 34.32 8.4 
10 2015 1007.7 402.5 395.6 29.87 34.23 6.5 
11 2015 1005.8 411.6 395.9 29.09 34.25 14.7 
12 2015 1006.9 400.0 399.0 30.52 34.12 0.9 
1 2016 1007.7 402.0 398.9 29.75 34.23 3.0 
2 2016 1007.2 409.5 398.4 29.76 34.34 10.5 
3 2016 1008.8 443.6 398.6 29.59 34.89 42.2 
4 2016 1008.7 476.6 399.8 28.85 35.12 73.6 
5 2016 1009.3 485.3 399.1 28.52 35.22 82.6 


Table S1. Monthly-averaged real time data from the TAO 0°, 170°W mooring used to calculate 
ApCO>. Monthly atmospheric pressure and seawater (sw) and air xCO> values are averaged from 
3-hourly MAPCO> data. Monthly SST and SSS values are averaged from daily-averaged 
National Data Buoy Center (NDBC) TAO CTD data (from: 
http://tao.ndbe.noaa.gov/tao/data_download/search_map.shtml). Note that the October 2014 
monthly average only includes data from 10/15/2014 - 10/31/2014; all other months include a 
full month of data. Surface CTD or real time data transmissions seemed to intermittently fail at 
this site between 11/16/2014 and 11/30/2015, and were replaced by data from the TAO 2°N, 
170°W mooring (during data overlap in the time series, TAO 2°N, 170°W SST and SSS were - 


0.04+0.44 (1o) and 0.16+0.18 (10), respectively, compared to 0°, 170°W). 
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